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Abstract 

Wave packet revivals and fractional revivals are striking quantum interference phenomena that 
can occur under suitable conditions in a system with a nonlinear spectrum. In the framework 
of a specific model (the propagation of an initially coherent wave packet in a Kerr- like medium), 
it is shown that distinctive signatures of these revivals and fractional revivals are displayed by 
the time evolution of the expectation values of physical observables and their powers, i.e., by 
experimentally measurable quantities. Moreover, different fractional revivals can be selectively 
identified by examining appropriate higher moments. 
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I. INTRODUCTION 



ution is gov- 



arising 



It is well known that the dynamics of a quantum wave packet whose evo 
erned by a nonlinear Hamiltonian can display a variety of non-classical effects 
from the interference between the stationary states comprising the wave packet. These ef- 
fects manifest themselves through several interesting observable phenomena such as revivals, 
fractional revivals and super- revivals of the wave packet at specific instants of timejsjj-Ql. 
Such revivals have been identified in several physical systems. The phenomenon has been 
observed experimentally and investigated theoretically in a wide class of suitably prepared 
states subject to specific nonlinear potentials, including Rydberg wave packets|3l|-|l0j| and 
molecular vibrational states^], 3|. For a recent comprehensive review of revival phenom- 
ena, see 

Detailed calculations have been made on revival times within the framework of the Jaynes- 
Cummings model of a two-level atom interacting with a single mode of electromagnetic 
radiation 



jl^. jl^]. Full and fractional revivals of angular momentum coherent states [r|. 
SU (1,1) coherent statesll7|, states evolving in specific one-dimensional potentials like the 

n n 

infinite square well [la], and general multi-level quantum systems[19] have also been studied. 

The non-classical effects that arise in the dynamics of specific wave packets evolving in a 
Kerr medium 2(|, collapse and revival phenomena displayed by a single-mode field propagat- 



ing in a nonlinear medium 



2l|, and the temporal behavior of a macroscopic superposition of 



coherent states governed by a nonlinear Hamiltonian |22j have been reported upon. An inter- 
esting analogy between revivals of wave packets and Poincare recurrences in classical maps 
has also been pointed out, and the interplay between revivals, recurrences and geometric 



phases has been elucidated 



24|. 



Much of the attention given to the revival phenomenon has been with a view to under- 
standing the precise nature of the transition from quantum to classical dynamics, and the 
deviations from classical predictions that are displayed during the long-time evolution of 
quantum states. Wave packets prepared so as to satisfy the minimum uncertainty condition 
at the initial instant of time are clearly appropriate candidates for this purpose. A Ryd- 
berg state (an electron excited to a high energy level in the bound-state spectrum by the 
application of a suitable laser pulse), for instance, is well described by a wave packet which 
is sharply peaked about a large quantum number uq, and the spectrum is approximately 
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linear in this energy region. Although such a wave packet generally spreads with time, it can 
display, under special circumstances, complete or full revivals as well as fractional revivals 
at specific instants of time during its evolution. 

What has emerged from both experimental and theoretical investigations of this phe- 
nomenon is the following: the broad features of wave packet dynamics are quite generic, 
regardless of the details of the physical system concerned, the initial wave packet consid- 
ered, and the specific nonlinear Hamiltonian governing its time evolution. Given an initial 
state \if}(0)), if the overlap function C{t) = \(ip(0)\ip(t))\ 2 returns to its initial value of unity 
at specific instants of time, we have a revival. Revival phenomena are essentially controlled 
by the parameters occurring in the first- and second-order terms in the Taylor expansion of 
the energy spectrum E n about the energy E no corresponding to the peak of the wave packet, 
namely, 



E — E 

J -'n ■ LJ no 

' - m\ 

m=l 



T — 

Z— J rn.\ 



d m E n 
dn m 

J n=no 



(n-n ) m . (1) 



The quadratic term in Eq. (0) can also lead to fractional revivals that occur at specific 
instants between two successive revivals of the wave packet. At these times the initial wave 
packet evolves to a state that can be described as a collection of spatially distributed sub- 
packets, each of which closely resembles the initial wave packet. Since each of the basis 
states comprising the wave packet acquires a different phase during its temporal evolution, 
in general, fractional revivals can occur only if the numerical values of the parameters that 
appear in the Hamiltonian satisfy certain conditions, namely, if certain ratios of these pa- 
rameters are sufficiently close to rational numbers. Anharmonic corrections due to the cubic 
and higher order contributions in Eq. JTJ, leading to so-called super- revivals, have also been 
.vested, in particular in the context of Rvdbe, ato mS B,fl,Q,0. For wave pac k - 
ets which are peaked sufficiently sharply about E no , these higher order effects that occur on 
comparatively longer time scales are negligible. 

While it is evident that revivals and fractional revivals are genuinely quantum mechan- 
ical in origin, it is of interest to examine the signatures of these phenomena in exper- 
imentally measurable quantities, i.e., upon the expectation values of appropriate physical 
observablesj3]. More specifically, we are concerned with the possibility of finding distinctive 
signatures of different fractional revivals occurring between successive revivals that would 
enable one to distinguish between them unambiguously. In this paper, we use the example 
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of a specific nonlinear Hamiltonian that is both physically relevant and for which analytic 
calculations can be performed, to show that this is indeed feasible. This implies that re- 
vival phenomena can be investigated directly, in terms of these experimentally accessible 
expectation values. 

The plan of the paper is as follows: In the next section we briefly outline the salient as- 
pects of revivals and fractional revivals relevant to our purposes, including the relationship 
between the time of occurrence of a specific fractional revival and the number of spatially 
distributed sub-packets representing the state at that instant of time. In Section III we 
obtain expressions for the mean and higher moments of the position and momentum opera- 
tors and track their behavior as the wave packet evolves in time. This helps us identify the 
distinguishing features displayed by these quantities at the instants when fractional revivals 
occur. We conclude with a discussion of our results in Section IV. 



II. REVIVALS AND FRACTIONAL REVIVALS 



Consider a nonlinear time-independent Hamiltonian H with a non-degenerate spectrum 
{E n } and eigenstates The initial state \ip(0)) in which the system is prepared is a 

superposition of these eigenstates, sharply peaked about some n . For ease of notation we 
shift n by no in Eq. and retain terms up to second order in r = n — no, since only these 
terms contribute to revivals and fractional revivals. The unitary time evolution operator 
U(t) = exp (—iHt/h) can be written in the natural basis in standard notation^] as 



•) , (2) 



apart from an overall phase factor. Here T c \ = h/\E'(no)\ and T rev = 2h/\E"(no)\. The 
quadratic term in the exponent is responsible for the spreading and distortion of the initial 
wave packet. However, if there exist specific instants of time such that, for every r in 
the summation, the exponent is equal to an integer multiple of 2iri, then U(t) reduces to 
the identity operator apart from an overall phase factor, and a revival occurs. Further, as 
mentioned in Sec. I, fractional revivals can occur at certain instants of time in between 
successive revivals. In general, exact revivals require special values of the ratio T c i/T rev , but 
generically C{t) can take on values arbitrarily close to unity at certain instants of time (2^. 
In some sense, such revivals may be regarded as analogs of Poincare recurrences of a classical 
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dynamical system in a coarse-grained phase space. 

For definiteness, we shall consider in the rest of this paper the example of wave packet 
evolution in a Kerr-like medium. The relevant physics of this situation is captured by the 
model Hamiltonian [2^,(2^] 

H = h X a ]2 a 2 = h X N(N - 1) (3) 

in the usual notation, with N = a) a, \ = positive constant. Essentially the same spec- 
trum models (^J the interaction energy of the atoms in a Bose-Einstein condensate. The 
initial state \ip{0)) is taken to be the oscillator coherent state \a) (where a G C), as it is 
a minimum position-momentum uncertainty state, and hence is "classical" in this sense. 
The corresponding position-space and momentum-space wave functions are Gaussian wave 
packets. Owing to the nonlinear nature of H, these wave packets not only spread out as 
t increases, but also change shape, and the state loses its coherence property. We have 
U(t) = Eo° ex P [~ix n { n ~ \ n ) ( n \) where \n) denotes the usual number operator eigen- 
state. Since n{n — 1) is always an even integer, it is evident that T d = T rev = 7r/x in this 
case. 

In between t — and t = T rev , fractional revivals occur at times t = 7tl/mx, where 
m = 2, 3, 4, . . ., and I = 1, 2, . . . , (m — 1) for a given value of m. This can be traced to 
the interesting periodicity property exhibited by U at these instants of time, consequent to 
which it can be Fourier expanded in the basis set {exp (— 2iril/m), < I < (m — 1)} for 
each m. Writing 

°° -i X n(n-l)t n 

\m)=e- ]al/2 J2- TT^^ W 



v n 

n=0 



we see that, at these instants of time, the exponential in the summand can be absorbed in 
the complex number a. Thus \ijj(7r/mx)) = iV'm) becomes a finite superposition of "rotated" 
coherent states with definite amplitudes, according to 

[ ET^fi^e- 2 ^), modd 

|V>m) = <^ (5) 

[ E^o 1 9i \ae i7r/m e ~ 2 * i[ l m ) , meven, 

where fi and gi are the Fourier coefficients. For instance, the state at time 7r/(2x) (corre- 
sponding to m = 2) is a superposition of the two coherent states \ia) and |— ia). In general, 
\iprn) is a superposition of m coherent states. The corresponding wave packet in position 
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space is a superposition of m spatially distributed Gaussian wave packets, and a fractional 
revival is said to have occurred. The periodicity property of U further implies that the wave 
function at times t = nl/mx, 1 < / < m — 1, is also a superposition of m wave packets (^J. 
It is therefore clear that experimental observations of fractional revivals that depend on 
probing the wave packet itself are quite involved and somewhat indirect, and careful inter- 
pretation of the results is required before firm conclusions can be drawn. It is this aspect 
that leads us to suggest that an examination of the moments of observables may provide a 
less intricate and more definitive means of identifying and sorting out revival phenomena. 



III. MOMENTS OF OBSERVABLES 



We now show how fractional revivals are mirrored, in distinctive ways, in the expec- 
tation values of the physical observables pertaining to the system. As the system enjoys 
revivals with a period T rev , all such expectation values are periodic functions of t with this 
fundamental period. 

The relevant observables are the position x = (a + aJ)/ a/2 and p — (a — a 1 )/iy/2. Clearly, 
their expectation values alone do not suffice to reproduce the full information contained in 
the wave function itself. In principle, an infinite set of moments, namely, the expectation 
values of all powers of x and p and their combinations, is required for this purpose. In this 
sense, the quantum system is equivalent to an infinite-dimensional classical system in which 
the role of the dynamical variables is played by the set of expectation values. However, we 
emphasize that even the first few moments can be seen to yield considerable information on 
the behavior of the system. 

Recalling that a\a) — a \a), it is convenient to define the c-number function 

a(t) = (^j(t) \a\i)(t)) = (a\ e lHt/h a e~ imlh \a) , (6) 

so that a(0) = a. In the special case in which H = H(N), this simplifies to the explicit 
expression 

oo I I 2r 

a(t) = a e-H 2 ]T ^- e -imr+i)-H(r)]t/H (7) 

n=0 r ' 

In general, this is a fairly complicated function of t. For the case at hand, however, with 
H = hxN(N — 1), some simplification occurs. We find 

a{t) = ae -H 2 (^s2 X t) [ cog Q a |2 sin ( 2x t)) _ z S i n (|«| 2 sin(2xt))] . (8) 
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Thus a(t) is a periodic function of time with period tt/x, as expected. 

It is convenient to introduce the notation a = a± + ia^ = {xq + ipo)/V% and v = \a\ 2 = 
\{xq +£>o). As is well known, x and p represent the locations of the centers of the initial 
Gaussian wave packets in position and momentum space, respectively. Then the expectation 
values of x and p can be written as explicit functions of t in the form 

(x(t)} = e-^a-co^x*) [ Xq cos(i/sin2x0 +Po sin (u sin 2 X 0] , ( 9 ) 



(p(t)} = e -Ki-cos2 X t) [_ Xo S i n ( us i n 2 x t)+p cos (u sin 2 X t)]. (10) 

Similarly, the second moments of x and p are found easily from the expression (a 2 (t)) = 
a 2 exp ( — 2ixt — v + z/e" 4 * x ') and its complex conjugate. We get 

2 (x 2 (t) ) = 1 + xl + pi + e"" (1 - cos 4xt) [(x 2 - pi) cos (2tf + v sin 4tf) 

+ 2rr p sin (2xt + ^ sin 4x0] > ( n ) 

2<p 2 (0> = 1 + xl + p\ - e-^- cos [(^-p 2 ) cos(2 X t + z/sin4 X 

+ 2rr p sin (2*t + ^ sin 4x0] • ( 12 ) 

For reasons already mentioned, the higher moments also carry much information of direct 
interest. The third moments can be written compactly in the form 

4<x 3 (0> = e-^- cos6 ^[(^-3^) cos(6x* + ^sin6xO + (3^po-piD 

x sin (6xt + v sin 6x0] + 6z/ [ M0) (1 + cos 2x0 + (p(0) sin 2 x0] , ( 13 ) 

4<p 3 (0> = e-^- cos6 ^[(^-3^) sin(6xt + ^sin6xO-(3^p -Po) 

x cos (6xt + v sin 6x0] + 6z/ [ (p(0) (1 + cos 2x0 - (^(0) sin 2 Xt] • (14) 

With these expressions and the explicit expressions for the fourth moments (which we do 
not write them down here), the variances of x and p as functions of t, as also the skewness 
and kurtosis in each case, can be obtained. The uncertainty product (or the product of 
the standard deviations) Ax Ap, which initially has the minimum value |, is of special 
interest. We do not write down the lengthy expressions for these quantities here, but we 
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shall comment upon their time variation in the sequel. Finally, we note that expressions for 
the higher moments can be deduced readily from the general result 

(a) k a k+l ) = a 1 i/Vf 1-008 2lxt) exp [-%x{l{l ~ 1) + 2kl) t - iv sin 2l X t] , (15) 

where k and I are non-negative integers. 

We now turn to a discussion of our results. This is most conveniently done with the help 
of plots drawn for typical values of the parameters concerned. 



IV. RESULTS AND DISCUSSION 



We first compare the explicit solutions found for the expectation values of x and p in 
Eqs. (JUJ) and (jl(Jj) with the solutions that would have been obtained for x(t) and p(t) had 
the system been a classical one, governed by the classical counterpart of the normal-ordered 
Hamiltonian H = hxa^ 2 a 2 , namely, H c \ = \{x 2 +p 2 ) 2 . Although the equations of motion 
corresponding to H c \ are nonlinear, it is evident that x 2 + p 2 is a constant of the motion, 
so that the phase trajectories are circles. However, the frequency of motion is dependent 
on the initial conditions (i.e., the amplitude of the motion), being equal to v — \[x\ +pl)- 
This is, of course, a well-known feature of nonlinear oscillators. But we note that the actual 
solutions for (x(t)) and (p(t)) in Eqs. © and (fTUj) are more complicated than the classical 
ones for x(t) and p(t) under H c \. This is a consequence of the quantum mechanical nature 
of the system, over and above the nonlinearity of H. However, the expressions for (x(t)) 
and (p(t)) can be given the following interesting interpretation in classical terms. Define the 
(non-canonical) pair of classical dynamical variables 

X = xe v{1 ' cos2xt \ P = pe u{1 - cos2xt) , (16) 

and the reparametrized time r = sin (2%t). The initial values Xq and Pq of these variables 
remain equal to Xo and po, respectively. Then Eqs. (JHJ and (JTOJ) can be re- written as 

X = Xq cos vt + Po sin vt , P = —Xq sin vt + Po c °s vr - (17) 

But these are the solutions to the system of equations 

dX „ dP 

^ = *P,- = -^, (18) 
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describing a nonlinear oscillator of frequency v = \(X% + P 2 ), in terms of the transformed 
variables (X, P) and the reparametrized time r. At the level of the first moments, therefore, 
the system is effectively a nonlinear oscillator after a suitable transformation of the relevant 
variables. 

Turning to the details of the time dependence of the various moments of x and p, there 
are two striking features that underlie the essential point we wish to make in this paper. 
First, the higher the order of the moment (or cumulant), the more rapid is its variation, 
since the leading frequency in the m th moment is 2m%. Second, the time dependence is 
strongly controlled by the factor exp \—v (1 — cos 2mxt)], m — 1, 2, ... , that modulates the 
oscillatory terms. While this permits substantial time variation for sufficiently small values 
of v, it acts as a strong damping factor for large values of u, except when cos (2mxt) is 
near unity. As one might expect, this happens precisely at revivals (when t = mr/x, an 
integer multiple of T rcv ), as expected. But it also happens, in the m th moment alone, at the 
fractional revival times t = (n + l/m)T TCV . Thus, by setting v at a suitably large value, we 
can ensure that the moments are essentially static, bursting into rapid variation at specific 
instants of time before reverting to quiescence. 

These points are illustrated by the figures that follow. Owing to an obvious symmetry of 
the Hamiltonian, the moments of x and p behave in an essentially similar manner, especially 
if we start with the symmetric initial condition x = p . Without loss of generality, we 
restrict ourselves to this case in what follows. We have set x — 5 in the numerical results 
to be presented, but this is irrelevant as all the plots correspond to t measured in units of 
n/X- We find that for very small values (<C 1) of xq and po (i.e., of u), the nonlinearity of 
H does not play a significant role, and the behavior of the system is much like that of a 
simple oscillator. Interesting behavior occurs for larger values of v. We therefore present 
results for three typical values of the parameters representing the initial conditions, namely: 
(a) x = po = 1 => v = 1; (b) x = p = VlO =>• v = 10; (c) x = p = 10 =>• v = 100. These 
correspond, respectively, to small, intermediate, and large values of v. In all the "phase 
plots" , the point representing the state at t — is labeled A. 

Figures l(a)-(c) show the variation of (x(t)) as a function of t for small, medium and 
large values of v. (As already mentioned, (p(t)) displays similar behavior.) For sufficiently 
large values of v, it is evident that, except for times close to integer multiples of T rev = n/x, 
(x(t)) and (p(t)) essentially remain static at the value zero. Figures 2(a)-(c) depict the 
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<x(t)> <x(t)> <x(t)> 

(a) (b) (c) 

FIG. 2: "Phase plot" of (p) versus (x) 

corresponding "phase plots" in the ((x) , (p)) plane. In Fig. 2(c), the representative point 
remains at the origin most of the time, except at times close to successive revivals, when it 
rapidly traverses the rest of the curve before returning to the origin. 

While sudden changes from nearly static values of (x(t)) and (p(t)) are thus signatures 
of revivals, the occurrence of fractional revivals is not captured in these mean values. The 
fractional revival occurring mid-way between successive revivals (e.g., at t = 7r/2x in the 
interval between t — and t = T rcv ), when the initial wave packet reconstitutes itself into 
two separate wave packets of a similar kind, leaves its signature upon the second moments. 
Figures 3(a)-(c) show the variation with time of the uncertainty product Ax Ap. In each 
case, this product returns to its initial, minimum value \ at every revival, rising to higher 
values in between revivals. Once again, for sufficiently large values of u, the product remains 
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FIG. 4: "Phase plot" of Ap versus Ax 

essentially static at the approximate value (| + v) for most of the time, but undergoes 
extremely rapid variation near revivals, and also near the fractional revivals occurring mid- 
way between revivals. During the latter, the uncertainty product drops to smaller values, 
but does not reach the minimum value |. 

There is a very striking difference in the behavior of the standard deviations near revivals 
as opposed to their behavior near the foregoing fractional revivals. This is brought out in 
Figs. 4(a)-(c), which is a "phase plot" of Ap versus Ax. For very small u, as in Fig. 4(a), 
Ax and Ap vary quite gently around a simple closed curve. When v is somewhat larger, as 
iri0Ja) which corresponds to v = 1, the plot begins to show interesting structure. For much 
larger values of v as in0Jc), the initial point A quickly moves out on the zig-zag path about 
the radial Ap = Ax line to the steady value represented by the point B, and returns to A 
at every revival along the complementary zig-zag path. Close to the fractional revival at 
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Xt/n Xt/ji Zt/ji 

(a) (b) (c) 



FIG. 5: Square of the skewness in x as a function of time 

t = (n+ |)T rcv , however, the representative point moves back and forth along the azimuthal 
path BCDB rather than the zig-zag path: clearly, a kind of "squeezing" occurs, as one of 
the variances reaches a small value while the other becomes large, and vice versa. (Of course 
the state of the system is far from a minimum uncertainty state throughout, except at the 
instants nT rcv .) 

The fractional revivals occurring at t = (n + |)T rcv and t — (n + |)T rcv , when the 
initial wave packet is reconstituted into a superposition of three separate wave packets, are 
detectable in the third moments of x and p. To make this unambiguous, we may consider 
the third moments about the mean values — or, in standard statistical notation, the square 
of the skewness, defined as 

P[ x) = ((x-(x)fy /({x*)-{xf)\ (19) 

and similarly for (3[ p \ Figures 5(a)-(c) show the variation of (3[ x ^ with t. It is evident that, for 

(x) 

sufficiently large values of v, (3{ remains nearly zero most of the time, except for bursts of 
rapid variation close to revivals and fractional revivals. Both p[ x ^ and (3^ actually vanish at 
t = nT TCV , but they remain non-zero at t — (n + //3)T rcv , I — 1, 2. More detailed information 
is obtained from a "phase plot" of (3^ versus (3[ x \ which we do not give here. 

Finally, we consider fractional revivals corresponding to m — 4, when four superposed 
wave packets appear. These are detectable in the behavior of the fourth moments of x and 
p. Equivalently, we may use the excess of kurtosis (/3 2 — 3), where the kurtosis of x is defined 
as 

^ = {(x-{x)Y)/{{x*)-{xf)\ (20) 
12 




Xt/jt Xt/jc Xt/jt 

(a) (b) (c) 

FIG. 6: The excess of kurtosis of x as a function of time 



with a similar definition for (5\ . The excess of kurtosis is the measure of the departure of 
a distribution from gaussianity. Figures Efa)-(c) depict how (f3 { 2 x) -3) varies with time. For 
sufficiently large z/, both (/^ — 3) and (/?2 — 3) remain essentially static near the value — | 
for most of the time. They vary rapidly near revivals, vanishing at t = nT rev because the 
wave packet is a Gaussian both in position space and in momentum space at these instants 
of time. As is clear from Fig. 4(c), they also vary rapidly near the fractional revivals at 
t — (n + Z/4)T rev (where / = 1, 2, 3), oscillating about the "steady value" — §. Once again, a 
"phase plot" of (/^ p) -3) versus (/^ — 3) (which we do not give here) helps identify features 
that distinguish between the three fractional revivals concerned. 

We have shown that distinctive, experimentally detectable signatures of the different 
fractional revivals of a suitably prepared initial wave packet are displayed in the expecta- 
tion values of physical observables and their powers. The complicated quantum interference 
effects that lead to fractional revivals can thus be captured in the dynamics of these expec- 
tation values, which may be regarded as the dynamical variables in a classical phase space. 
While this is, in principle, an infinite-dimensional space, what is relevant in practice is the 
temporal behavior of a finite number of moments, since fractional revivals corresponding to 
very large values of m are not easy to detect in any case. 
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